Josephson Effect between Condensates with Different Internal Structures 
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A general formula for Josephson current in a wide class of hybrid junctions between different inter- 
nal structures is derived on the basis of the Andreev picture. The formula extends existing formulae 
and also enables us to analyze novel _B-phase/A-phase/B-phase (BAB) junctions in superfluid 3 He 
systems, which are accessible to experiments. It is predicted that BAB junctions will exhibit two 
types of current-phase relations associated with different internal symmetries. A "pseudo-magnetic 
interface effect" inherent in the system is also revealed. 
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Superfluid systems with internal degrees of freedom 
produce diverse ordered structures, which provide a new 
arena for exploring the fertile physics behind them. The 
Josephson effect [Q extracts the global phase resulting 
from the spontaneous breakdown of gauge symmetry. 
Moreover, it also provides us with information on the 
internal structures formed from other broken symmetries 
since quasiparticles acquire an excess phase from the in- 
ternal structures while traveling through the junction. In 
fact, so-called n junctions associated with high-T c super- 
conductors provide convincing evidence for d-wave sym- 
metry and recent experiments on Pb-S^RuGi- 
Pb Josephson junctions M could be explained by p-wave 
symmetry in Sr 2 Ru04 [|| o. Moreover, metastable it 
states observed in superfluid 3 He weak links |7| [| |) are 
considered to be a signature of texture due to the internal 
degrees of freedom of p-wave order parameters jl0|, [ll] . 
In the above ways, order parameter configurations have 
been investigated by using Josephson current-phase rela- 
tions. 

Since a wide variety of structures could be realized 
in systems with internal degrees of freedom such as the 
textures in superfluid 3 He systems, we require a general 
theory for the Josephson effect to enable us to identify 
the structures. In this Letter, we derive a general for- 
mula for Josephson current, applicable to a wide class of 
hybrid junctions between different internal structures in 
unitary states, on the basis of the Andreev picture p2l of 
the Josephson effect developed by Furusaki and Tsukada 
(FT) [ pUf , including explicit expressions for Andreev re- 
flection coefficients reflecting their structures. 

Consider a junction composed of three different super- 
fluid regions separated with flat interfaces perpendicular 
to the z-axis. The interfaces are located at z = and 
z = L. We shall denote these three regions by L (z < 0), 
C (0 < z < L), and R (L < z), respectively. We assume 
that the order parameter is uniform in each region. We 
also assume that the effective mass m, Fermi velocity vp, 
and Fermi wave number Uf are the same in all regions. 
The potential barriers at the interfaces are ignored for 



simplicity. We shall use a triad (x,y,z) as the base for 
a frame of reference. 

The dc Josephson current I can be calculated from 
the temperature Green's function, Q(z, z', kn,u> n ), which 
is obtained by an analytic continuation, E + i0 + — > iu> n , 
from the retarded Green's function constructed by scat- 
tering data. Here u) n is the Matsubara frequency defined 
by ut n — 7r(2n+l)//3 (n = 0, ±1, ±2, ...) at temper- 
ature T = and feii = (k x ,k y ,0) is defined from 
k = (k Xl k yi k z ) (\k\ = 1) along the classical trajectory. 
Note that the Green's functions are extended in the form 
of a 4 x 4 matrix due to the spin and particle-hole de- 
grees of freedom required for triplet pairing. Taking the 
particle number conservation into account, the current at 
z = is calculated by 

Let us calculate the scattering data necessary for con- 
structing the temperature Green's function. In a quasi- 
classical approximation, the motion of a quasiparticle of 
energy E is described by Andreev equations applicable 
to an arbitrary type of pairing fll4| as 

((-\hv F k z §- z -E)a A* \/ \ =() 

^ A^t (UtoFk,&-E)v ) W ' 

(2) 

where u = (lit u^f and v — (u-f fj.)'. <J is the 2x2 
unit matrix. The gap matrix in spin space is denoted by 
A k . In each region (a = L,C,R) the gap matrix, A£, is 
assumed to be constant and unitary: A^-A^ = |A^| 2 <7 . 
The solutions of the Andreev equations in each region are 
expressed as 
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= M*{z)S, (3) 
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where fi£ = yj E 2 - |A*| 2 , S p = (S pT S pl )\ and S h = 

(Sh-\ ShiY- S p (h)cr with cr =t or J, is the probabil- 
ity amplitude for the mode propagating as a particle-like 
(hole-like) quasiparticle with spin cr. The two amplitudes 
<Sl and <Sr are connected by the continuity equation: 

5 L = Mf(0)- 1 M5(0)M5(L)- 1 M|(i)5 R . (4) 

Here we introduce scattering coefficient matrices in the 
following form to express the scattering processes con- 
cisely; 



S 1T s n 



(•5) 



where s aa i is the scattering coefficient for the process in 
which an injected quasiparticle with spin a' is scattered 
as a quasiparticle with spin cr. The retarded Green's 



function is constructed from the scattering coefficient ma- 
trices under the FT prescriptions |u|. For z < z' < 0, 
the retarded Green's function is given by 

G{z,z',k\\,E) = -i 



h 2 kpk z 



(E + n 



k- i 



(6) 



where is the 2x2 zero matrix and fc_ = ku — k z z. The 
Andreev reflection coefficients for a process in which a 
particle-like quasiparticle injected from L is reflected as 
a hole-like quasiparticle is obtained by solving Eq. M) as 
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where f* = A*/(E + fi*) and ^ c = tl%/hv F k z . The 
matrix &2 describes the reverse process. The matrices 
0,2 and fii are related by a^iE, k) = a\[E, fe_)'. Substi- 



tuting the temperature Green's function obtained by the 
analytic continuation into Eq. ([!]), a general formula for 
the Josephson current at temperature T is derived as 



I =- 



d 2 h 
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where f2^ L = y ^ + | A£ | 2 • This is the central result 
of this Letter. This formula is applicable to any type of 
Josephson junction between unitary states with any sym- 
metry. Note that our formula still preserves the original 
FT form expressed by the difference between Andreev re- 
flection coefficients, describing the net current carried by 
the two processes: the scattering of a particle-like quasi- 
particle into a hole-like quasiparticle and its reverse pro- 
cess. This formula can indeed be reduced to the FT for- 
mula by substituting the s-wave order parameters. Our 
formula is, therefore, a natural extension of the FT for- 
mula. 

Let us consider certain special cases reported previ- 
ously by several authors. Substituting the order param- 
eters used in Ref. [TTj into our formulae (^) and (||), and 



adopting the limit L — > 0, we obtained an analytic ex- 
pression for the Josephson current through a pinhole be- 
tween two 3 He-B reservoirs (Eq. (2) of Ref. [k]). This 
reduces to the Kurkijarvi formula |H| for the B phase if 
the n vectors of two 5-phase reservoirs are parallel. This 
is the same as the expression for the s-wave supercurrent 
through a short orifice fig] . By employing the ^4-phase 
rather than the 5-phase order parameter, the Kurkijarvi 
formula for the A phase is reproduced. Let us take 
another limit Ac — > with a finite L. This reproduces 
the current-phase relation of the superconductor/normal- 
metal/superconductor (SNS) junctions without a poten- 
tial barrier ]Tfl] . The Tanaka-Kashiwaya formula 
for unconventional singlet superconductors without a 
potential barrier is also reproduced. In addition, the 
Yamashiro- Tanaka-Kashiwaya results || are reproduced 
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when we employ their order parameters. Thus our for- 
mula covers the previous formulae for concrete examples. 

Furthermore, our formula enables us to analyze new 
kinds of junction. One example is B-phase/ ^4-phase/i?- 
phase (BAB) junctions in superfluid 3 He systems. This 
type of junction could be realized by using the experi- 
mental setup used to study the Andreev reflection E^] 
at AB interfaces. It is also suggested that a pseudo- 
A phase could be formed around the orifice in 3 He-i? 
weak links pc[ |. In spin-triplet Cooper-pair condensates, 
the boundary condition should be carefully considered 
to include spin-orbit degrees of freedom. The order pa- 
rameter of the A phase is defined by a triad (w±, W2, 1) 
in orbital space and a vector d in spin space. For the 
B phase we need a rotational matrix R(n, 9) with ro- 
tational axis n and rotational angle 9 relating the spin 
space to the orbital space. We take the phase angle for 
the A phase to be zero. For an AB interface with dipole 
energy |2l], ^2|, the boundary conditions are 1 || d J_ z, 
d, = R(n,9L)wi, and i&i || z, where 6l — cos _1 (— 1/4) 
(Leggett angle). Since the A phase always has the am- 
biguity of d — > — d, Wi — > —Wi, and w>2 — > —102, we 
choose W\ — z, 11)2 — x, and I — y without loss of gener- 
ality. The vector n in the B phase is fixed in one of four 

directions, i.e., 1—^/3/5, ±^/l/5,±^/l/5J for d = I and 



^ v /375 j Tn /I75, ±yfTj5 ) for d = I. We use the letters 

a, b, c, and d in turn to specify the different orientations 
of n for brevity fic|| . The order pairs of the letters aa, 
ab, etc., thus denote the order parameter configurations 
at the two interfaces. The order parameters for the con- 
figurations of aa and ab, which become important in the 
following discussions, are given by 

A£ = A B e- [ip/2 /2 

x |2i(T fc 2 + (7i(VSk x - k y ) - &z(k x + n/SA^) j 

A£ = A A (k x - ik z )a 
A| = A s e^ /2 /2 
x |2i(T fc z + (7i (±V3k x - k y ) - a 3 (k x ± V3ky)} , (9) 

where ip is the phase difference and (<7i, 02, 03) are Pauli's 
matrices. A a and Ab are the gaps for the A phase and 
the B phase, respectively. A rotation of ir around the y- 
axis affects the transformations a<->d, b^>c, L<-»R, and 
d — * — d, where we used the ambiguity of the A phase. 
In addition, aa and bb configurations give the identi- 
cal current, same as ab and ba. As a result, we pre- 
dict two types of current-phase characteristics: type (a) 
{aa,bb,cc,dd} and type (b) {ab,ba,cd,dc}. Figure |l| 
shows the temperature dependence of the Josephson cur- 
rent in (a) aa and (b) ab configurations for L = 2£o 
where the coherence length £0 in the B phase at zero 
temperature is given by £0 = [7C(3)/48] 1/2 (^ f /7tA:bT c ) 
with ( being the Riemann zeta-function. Here we assume 




12 3 4 

Phase Difference <p/% 



12 3 4 

Phase Difference ip/rc 



FIG. 1: Josephson current-phase relations of BAB junctions 
for L = 2£ at temperatures 0.1T C , 0.3T C , 0.5T C , 0.7T C , and 
0.9T C in decreasing order of gradient at the origin for the con- 
figurations of (a) aa and (b) ab. Here, To = 5 , As(0)fc|'/47rft 
with S being the interface area. 



4* 0.5 





12 3 4 

Phase Difference <p/tt 



12 3 4 

Phase Difference ip/rc 



FIG. 2: Josephson current-phase relations of BAB junctions 
at O.lTc for L = 0, £0, 2£o, 3£o, 5£o, and 1000£o, in decreasing 
order of gradient at the origin for the configurations of (a) aa 
and (b) ab. 



the gaps obey the BCS temperature dependence. We can 
observe current-phase curves with 47r rather than 2-7r pe- 
riodicity and find some peak structures at lower temper- 
atures. In particular, the structures become pronounced 
as the temperature decreases. In the Andreev picture, 
the Josephson current is carried by both a discrete part 
(Andreev bound states) and a continuous part of the exci- 
tation spectrum. This suggests that the peak structures 
and the 4-7T periodicity come from the Andreev bound 
states and the continuous spectrum, respectively. 

It is known that Andreev bound states are localized 
near an interface. In the system with two interfaces that 
we consider here, the interaction between bound states 
near the two interfaces results in the current depending 
on the distance between the interfaces |2^, [m). In order 
to clarify the bound state contribution to the current, we 
investigate the current-phase relations of junctions with 
different ^4-phase widths at low temperatures. Figure ^ 
clearly shows that the peak structures are pronounced, 
while the 47r component is suppressed as L decreases. 
When L ^> £ > the Andreev bound states near the two 
interfaces are isolated and their overlap is lost. A BAB 
junction is regarded as a series of BA and AB junctions. 
Since the total phase difference through the BAB junc- 
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tion is (p, the BA (AB) junction has ip/2 dependence, 
resulting the 4-7T periodicity in the current-phase charac- 
teristics. In addition, the current-phase relations in (a) 
and (b) merge into a single curve in the limit L — > oo be- 
cause the current contributions from isolated interfaces 
are the same. When L — > 0, our model is reduced to the 
pinhole model in 3 Hc-B under a sufficiently large mag- 
netic field applied parallel to the wall |l0| . The surface 
free energy — (if • Rz) 2 produced by the applied mag- 
netic field H is minimized if fi aligns in one of the four 
directions, a, b, c, or d. Therefore, I plays a similar 
role to H in determining ft directions. This could be re- 
garded as being a "pseudo-magnetic interface effect" that 
is inherent in BAB junction systems. The difference be- 
tween H and the AB interface appears in the number 
of allowed h directions: all four directions for H, while 
only two directions in each d = ±1 for the AB interfaces. 
Thus the curve for L = in Fig. ||(b) is simply Yip's 
7r state Jl0| , that is produced by the current cancellation 
due to quasiparticles with different n-dependent effective 
phases. The curve in Fig. 0(a) corresponds to the curve 
for "AA" in Fig. 4 of Ref. [l(J. From the above discussion, 
the curve for "AC" cannot exist in a BAB junction. For 
L ~ £cb the Andreev bound states near two interfaces 
maintain their overlap and the BAB junction behaves 
like a single junction. The current-phase relation still 
maintains the characteristics of the L — curve. 

Now let us revisit Fig. |l|. The bound states localized 
near each interface overlap because L ~ £o- Since the 
current contribution due to the Andreev bound states 
is dominant at low temperatures, 27r-periodic current- 
phase characteristics appear in the curves as peak struc- 
tures. As the temperature becomes higher, the contin- 
uous spectrum contribution becomes large and obscures 
the Andreev bound state contribution. Therefore, we 
conclude that the peak structures appearing in Fig. [l] 
come from the Andreev bound states and the 4-7T period- 
icity is brought about by the continuous spectrum. 

Throughout the above studies of BAB junctions, the 
current-phase characteristics reflected the configurations 
of textures and were modified by arranging the textures 
of each phase. Since our formula applies to any type of 
Josephson junction, systematic studies of current-phase 
relations reveal the internal structures due to spin-orbit 
degrees of freedom as well as their pairing symmetries. 

In summary, we derived a general formula for Joseph- 
son current based on the Andreev picture. This formula 
made it possible to analyze any type of Josephson junc- 
tion between unitary states with any symmetry, and was 
applied to novel BAB junctions in superfluid 3 He, which 
are accessible to experiments. It has been predicted that 
there are two types of current-phase relations associated 
with AB interfaces inherent in the system, acting sim- 



ilarly to magnetic field H through I in the A phase 
("pseudo-magnetic interface effect"). 
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Note added. — During the preparation of this Letter, 
we learned of a theory |25| related to Josephson current 
in a similar system. 
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